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NUMERICAL  MODELING  OP  AIRBLAST 


CONTRACT  No.  N00014 -86-C-2197 


1st  Year  Final  Report 


Introduction 

The  SAIC  effort  during  the  period  March  10,  1986  to 

March  10,  1987  is  described.  Included  in  this  report  are 
appendices  for  the  deliverables  such  as  fortran  code 
listings.  The  effort  was  scaled  back  due  to  budget 

considerations.  The  areas  that  are  contained  herein  are 
Tasks  2.1,  2.2,  2.4,  2.5,  2.6,  2.7,  2.8,  and  2.9  as 
described  in  the  statement  of  work.  SAIC  over  the  last 
year  has  provided  a  highly  qualified  staff  of  research 
scientists  to  perform  the  tasks  described  under  the  contract 
statement  of  work.  Dr.  Mark  A.  Fry  has  acted  as  principal 
scientist  for  the  contract  work.  Dr.  Shmuel  Eidelman  has 
worked  on  site  for  8  months  and  then  4  months  offsite  at 
SAIC.  Mr.  Ken  Laskey  assisted  as  a  junior  level  research 
scientist  for  3  months.  Mr.  Mike  Zuniga  joined  the  effort 
in  November  1986  as  a  junior  level  research  scientist. 

1 
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The  2-D  axisymetric  FCT  code  was  applied  to  nuclear  and 
non-nuclear  problems.  A  burn  routine  for  HE  has  been  added 
to  the  code  but  has  not  been  thoroughly  checked.  The 
Minor-scale  event  was  simulated  and  waveforms  at  various 
station  locations  were  obtained.  These  were  used  in  a 
comparison  with  a  full  3-0  simulation  and  with  the 
experimental  data  it  was  determined  that  the  simulated 
thermal  layer  produced  multiple  peaks  in  the  overpressure 
waveforms.  This  was  primarily  due  to  a  tear  in  the  plastic 
bag  used  to  contain  the  helium.  The  numerical  simulations 
did,  however,  agree  quite  well.  The  results  of  this  work 
were  presented  at  the  Minor  Scale  Symposium  in 
Albuquerque . 

The  2-D  code  was  further  optimized  for  running  on  the 
CRAY  class  computers  by  employing  pointers  in  the  data 
structure.  Speed  increases  of  10  to  15%  were  obtained. 

Work  was  also  performed  in  addressing  improvements  for 
data  management  in  2  -D  codes.  This  work  has  been  titled 
"Monotonic  Logical  Grid  Algorithms”.  Mr.  Mike  Zuniga  has 
been  working  full  time  on  these  improvements. 


A  simulation  program  has  been  written  to  test  the  algorithm. 
In  addition*  a  graphics  program  has  also  been  written  to  aid 
in  this  research. 


TASK  2.2 

A  new  first  principles  model  for  dust  ingestion  was 
developed.  The  momentum  integral  equation  is  solved  in  ID 
along  the  ground  surface.  This  equation  along  with  boundary 
layer  physical  assumptions  provides  an  unsteady  prediction 
of  the  mass  flux  emerging  from  the  ground.  This  model  was 
first  proposed  by  Mirels  and  is  now  implemented  into  the  2- D 
FCT  code.  The  option  D(JST=2  initiates  the  use  of  this 
model.  A  listing  of  the  2 -D  code  and  the  subroutine*  DUSTN, 
can  be  found  in  Appendix  A. 

TASK  2.4 

The  3-D  PCT  code  was  used  to  study  the  late  time  cloud 
rise  geometry  from  multiple  nulear  explosions.  Calculations 
involving  two  bursts  separated  by  distance  and  time  were 
executed.  Cases  of  1  MT  bursts  at  20  and  30  second 
separation  with  distance  separation  of  250  m*  500  m*  1  km. 
and  3  km  were  completed.  The  results  indicate  a  strong 
sensitivity  in  cloud  formation  with  a  small  distance 


separation.  Previous  system  modeling  codes  assumed  no 
sensitivity  and  have  now  been  corrected. 

Spike  attacks  involving  many  bursts  simultaneously  in  a 
20  x  5  km  grid  were  investigated.  Simulations  performed 
with  the  3-0  code  indicated  a  strong  asymmetry  out  to 
several  minutes  in  time  and  up  to  40  to  50  km  in  altitude. 
A  paper  summarizing  the  results  was  presented  at  the  DNA 
DUST  Symposium  in  July  ,  1986.  The  abstract  of  the  paper  is 
presented  in  Appendix  B.  Additional  presentations  were  made 
at  DNA  headquarters. 

As  an  additional  item,  particles  were  added  to  the  3  -D 
code  for  tracking  material  with  time.  The  3-D  code  was  used 
to  study  the  possiblity  of  using  distributed  chemical 
charges  in  array  geometries.  Numerous  chemical  energy 
release  models  where  employed.  The  results  are  summarized 
in  a  report  in  Appendix  C. 

TASK  2.5 


Work  on  this  task  was  not  performed  due  to  budget 
limitations. 
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TASK  2.6 


Work  on  this  task  was  not  performed  due  to  budget 
limitations. 

TASK  2.7 


Work  on  this  task  was  not  performed  due  to  budget 
limitations. 

TASK  2.8 

The  3-D  graphics  program  has  been  updated  and  improved 
for  hard  copy  display  on  Tektronix  compatible  terminals. 
Links  were  written  to  use  DXSSPLA  as  well  as  the  NCAR 
graphics  package.  A  choice  of  bold  face  lettering  is  now 
also  available.  The  extensive  changes  and  the  improved 
documentation  are  available  in  the  source  listing.  Appendix 
D  . 

TASK  2.9 


Task  2.9  was  not  performed  due  to  budget  considerations. 


APPENDIX  A 


SOURCE  LISTING  OP  THE  2-D  AXISYMMETRIC 
AND  CARTESIAN  NUCLEAR  AND  NON-NUCLEAR 


MODELLING  CODE 
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Three-Dimensional  Effects 
in  Dust  Clouds-  A  Numerical  Approach 

Mark  A.  Fry 

Science  Applications  International 

Three-dimensional  numerical  simulation  codes  combining  minimal  diffusion  algorithms 
which  are  exceptionally  fast  on  today’s  supercomputers  have  been  used  to  investigate  nu¬ 
clear  dust  cloud  scenarios.  Results  show  persistent  multi-dimensional  effects  for  long  peri¬ 
ods  of  time.  Conclusions  based  solely  on  two-  dimensional (2-D)  calculations  are  misleading 
and  sometimes  incorrect.  Modeling  of  high  explosive  dust  cloud  scenarios  has  also  been 
performed  and  a  unique  high  explosive  and  numerical  simulation  program  for  dust  clouds  is 
proposed.  When  more  than  one  explosion  occurs  with  proximity,  multi-  dimensional  effects 
become  important.  Separation  distances  up  to  many  fireball  radii  produce  highly  complex 
flows  which  result  in  asymmetrical  dust  cloud  geometries.  Even  separation  distances  less 
than  a  fireball  radius  can  influence  cloud  shape.  The  need  to  perform  parametric  analysis 
of  these  effects  has  led  to  the  development  of  a  very  feist  non-diffusive  3-D  code,  FAST3DD 
with  dust  entrainment.  Pertinent  aspects  of  the  physics  modeled  will  be  explained.  High 
explosive  dust  cloud  simulation  provides  a  research  path  when  combined  with  numerical 
calculations  can  lead  to  new  understanding  of  entrainment  processes  and  cloud  develo- 
ment.  Improvements  in  the  understanding  of  the  dust  lofting  mechanism  have  been  made 
with  implementation  of  a  turbulent  of  a  turbulent  boundary  layer  model  for  use  with  in 
2  and  3-D  codes.  Multi-cloud  experiments  with  HE  are  proposed  with  direction  from  3-D 
calculations.  Spherical  charges  of  1000  pounds  or  more  that  can  be  remotely  detonated 
and  launched  into  a  proper  configuration  would  be  used.  A  sample  scenario  has  been 
computed  and  results  will  be  shown. 
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Introduction 


This  report  describes  the  results  of  the  one  year  feasibility  study  of  the  Dis¬ 
tributed  Explosive  Charge  Array  (DECA)  concept,  as  applied  to  antitank  mine 
neutralization.  The  DECA  concept  was  proposed  as  an  area  weapon  which  would 
have  the  maximum  blast  effect  for  a  given  weight  of  explosive.  The  goal  of  our  study 
was  to  substantiate  this  claim,  by  detailed  numerical  simulations  and  blast  exper¬ 
iments.  In  the  current  study  we  considered  DECA  arrays  formed  by  line  charges 
of  explosive  positioned  with  selected  distances  between  charges.  The  simulation  of 
the  blast  wave  propagation  from  a  multitude  of  finite  length  explosive  charges  was 
done  on  a  CRAY-XMP  supercomputer,  using  a  fully  3-Dimensional,  time  depen¬ 
dent  numerical  model.  Parametrical  study  of  the  maximum  pressure  and  impulse 
dependence  on  surface  density  was  done  for  DECA  formed  from  strands  of  Prima- 
cord.  The  formation  and  propagation  of  blast  waves  from  the  DECA  was  studied 
for  arrays  suspended  above  the  ground,  and  lying  on  the  ground,  with  different 
spacings  of  the  Iremite  or  Primacord  charges. 

Experiments  were  conducted  at  SRI  International  remote  test  facility  at  Corral 
Hollow,  California.  Experiments  were  done  with  full  size  DECA  arrays  of  Iremite 
or  Primacord  explosives,  and  were  designed  to  follow  numerical  simulations.  The 
►  experimental  data  was  used  to  verify  the  computer  simulations.  To  illustrate  the 

effectivness  of  DECA  against  antitank  mines,  in  two  tests  M-15  mines  were  desroyed 
by  DECA  arrays,  suspended  25cm  and  50cm  above  the  ground. 
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1.  Numerical  Simulations 

For  numerical  simulation  of  the  blast  waves  of  complex  structures  generated  by 
DECA  explosions,  we  used  FAST3D  computer  code.  This  code  was  developed  re¬ 
cently  in  the  Laboratory  for  Computational  Physics  at  NRL,  and  it  solves  unsteady 
and  fully  three-dimensional  flow  problems,  with  good  resolution  of  the  shock  waves 
and  contact  discontinuities.  This  code  is  vectorized  and  optimized  for  the  most 
efficient  use  on  the  Cray  X-MP  supercomputer.  More  details  about  this  code  struc¬ 
ture  and  the  numerical  method  used  are  given  in  Ref.l.  This  reference  is  added  as 
Appendix  A  to  the  current  report. 

The  purpose  of  the  current  study  was  to  achieve  levels  of  overpressure  and 
impulse  loads  which  cause  collapse  of  the  antitank  mine  casing.  According  to  the 
BRL  study  (Ref  2.),  published  in  19S3,  pressure  levels  of  above  2000psi  and  the 
blast  wave  impulse  of  above  950psi  •  sec  will  cause  rupture  of  the  M-15  mine  casing, 
in  the  area  of  the  central  fuse  cavity.  This  datum  was  taken  as  a  reference  point  for 
defining  the  range  of  explosive  surface  density  able  to  produce  the  load  needed  to 
destroy  antitank  (AT)  mines. 

The  size  of  the  DECA  array  was  first  choosen  to  be  4m  x  4m.  This  is  a  relevant 
size  for  potential  applications,  and  it  allowed  us  to  estimate  the  magnitude  of  the 
edge  effects  (nonplanar  relief  of  the  blast  energy  at  the  edges  of  the  DECA  charge). 
Later,  it  was  realized  that  the  edge  effects  are  small  at  the  time  when  most  of 
the  blast  load  affects  the  ground,  and  it  was  decided  to  study  arrays  of  2m  x  2m  in 
order  to  reduce  expenses  of  the  experiments,  and  improve  the  accuracy  of  numerical 
simulations. 

In  order  to  assure  good  resolution  of  the  blast  waves  in  all  spatial  directions,  all 
the  simulations  were  done  on  an  evenly  spaced  computational  mesh  (60  x  50  x  50). 
The  linear  geometry  of  the  individual  charges  was  chosen  because  of  the  simplicity 
of  its  implementation,  both  in  experiments  and  simulations.  Only  one  fourth  of 
the  physical  plan  was  simulated  numerically,  because  of  symmetry  considerations. 


For  a  typical  simulation,  at  the  time  t  =  0  we  deposited  energy  and  densitv  i  cor¬ 
responding  to  levels  observed  for  the  explosives  under  consideration)  into  cells  of 
the  computational  domain,  located  along  straight  lines  and  evenly  spaced.  Then 
an  inital  value  problem  was  solved  explicitly  for  the  system  of  three-dimensional 
Euler  equations.  The  pressure,  density,  velocities  (in  X,Y  and  Z  directions)  and 
energy  were  recorded  in  the  form  of  binary  arrays  of  values  for  every  grid  point 
in  the  computational  domain.  Because  of  the  large  size  of  these  arrays,  they  were 
recorded  only  for  a  given  increment  of  time  or  integration  step  (usualv  S  times  per 
simulation).  However,  all  physical  parameters  were  recorded  at  every  integration 
step  for  144  points,  located  in  the  most  relevant  places  of  the  computational  domain 
(including  the  points  corresponding  to  the  location  of  the  pressure  transducers  in 
the  experiments). 

For  analysis  of  the  complex  three-dimensional  flow  field  developing  after  a  typ¬ 
ical  DECA  explosion,  a  package  of  graphic  subroutines  was  developed.  The  plots 
produced  by  these  subroutines  are  demonstrated  below,  where  we  discuss  the  results 
of  our  simulations.  Here  we  list  the  various  functions  of  the  graphics  routines: 

1.  Contour  plots  of  X-Y  cross  sections  of  the  domain; 

2.  Contour  plots  of  X-Z  cross  sections  of  the  domain; 

3.  4.  The  same  as  1.  &  2.  in  color; 

5.  Time  history  of  pressure,  impulse,  density  and  temperature  (or  total 

pressure); 

6.  Contour  plots  of  .the  impulse  distribution  on  the  ground; 

7.  Graphs  of  pressure  and  impulse  distribution  on  the  ground. 

The  first  set  of  computations  and  experiments  was  done  with  Iremite  explo¬ 
sive.  The  surface  density  of  the  explosive  was  about  2.6kg/m2.  Figure  1.1  shows 
a  schematic  layout  of  a  typical  simulation.  Initially,  the  explosive  arrays  consisted 
of  nine  4m  long  lines  of  Primacord  spaced  0.5m  and  placed  on  the  ground.  Be¬ 
cause  of  symmetry,  in  the  computational  domain  we  have  4.5  lines  (the  line  at  the 
center  of  the  array  is  cut  in  half  by  the  symmetry  plane  Z-Y)  and  the  lines  are 


2m  long  (again,  physical  lines  of  explosive  are  cut  in  half  by  the  symmetry  plane 
Z-X  ).  In  Figures  1.2(a,b,c,d),  pressure  contours  are  shown  in  the  X-Z  plane  for 
the  blast  wave  propagating  after  the  explosion  of  a  1001b  charge  of  DECA  explosive 
formed  by  nine  4m  charges  of  Primacord.  The  pressure  contours  are  shown  for 
t  =  0.13ms,  t  =  0.30ms,  t  =  0.40ms,  t  =  0.78ms  after  the  detonation.  It  is  notice¬ 
able  that  in  this  case,  because  of  the  0.5m  space  between  the  charges,  the  load  on  the 
ground  is  not  uniform.  Also,  the  edge  effects  are  significant  in  this  case.  In  Figures 
1.3(a,b,c,d),  pressure  contours  are  shown  for  the  X-Y  plane  (the  ground).  The  max¬ 
imum  pressure  in  these  plots  is  allocated  along  the  lines;  although  this  maximum 
is  moving  along  the  ground  in  time  (compare  Fig.  1.3b  and  Fig.  1.3c),  this  is  not 
enough  to  issure  an  even  load  on  the  ground.  In  Figures  1.4(a,b,c)  the  contour  plots 
of  the  impulse  on  the  ground  are  shown  for  t  =  0.25ms,  t  =  0.50ms,  t  =  1.0ms,  after 
detonation.  The  maximum  impulse  achieved  in  this  simulation  was  %  ISOOpsi  •  ms. 
These  levels  were  observed  directly  under  the  lines  of  explosives.  In  the  areas  be¬ 
tween  the  lines  of  Primacord,  the  impulse  is  around  450psi  •  ms.  These  areas  of 
small  impulse  are  s:  20cm  wide.  The  blast  wave  contibutes  most  significantly  to  the 
impulse  in  the.  first  0.5ms  after  the  detonation. 

In  Figures  1.5(a,b),  graphs  of  pressure,  impulse,  temperature  and  density  are 
shown  for  the  same  simulation,  as  functions  of  time  elapsed  after  the  detonation. 
(The  coordinates  of  the  station  where  the  recording  was  made  are  given  in  the  third 
line  of  the  graph  headings.)  The  station  coresponding  to  Fig.  1.5a  is  located  under 
the  explosive.  That  leads  to  very  high  presssure  and  impulse.  However,  for  the 
station  located  between  the  lines  of  explosives  shown  in  Fig.  1.5b,  the  maximum 
pressure  is  only  about  l.Skpsi  and  maximum  impulse  is  500psi  •  ms.  Initially,  the 
uneven  load  distribution  was  not  considered  a  negative  factor,  since  it  was  assumed 
that  to  destroy  an  AT  mine  it  would  be  sufficient  to  produce  a  very  high  load  on 
part  of  the  mine  .  Later,  it  was  learned  that  this  concept  could  be  applied  only  to 
AT  mines  which  are  activated  by  large  pressure  plates,  and  it  would  be  inefficient 
for  blast  resistant  mines.  For  this  reason,  later  in  our  study  we  tried  to  get  more 


even  distribution  of  the  blast  wave  load  on  the  ground. 

Let’s  consider  in  more  detail  the  dynamics  of  the  blast  wave  for  DECA.  Simul¬ 
taneous  detonation  of  a  Distributed  Explosive  Charge  Array  formed  by  line  charges 
of  explosive  will  at  first  lead  to  cylindrical  blast  waves  in  the  vicinity  of  the  charges. 
At  a  later  time,  cylindrical  blast  waves  from  single  charges  of  the  array  collide  and 
after  multiple  reflections,  a  pseodo  planar  wave  would  form.  Since  the  main  advan¬ 
tage  of  the  DECA  concept  is  the  use  of  planar  blast  waves  to  efficiently  produce 
blast  loads  on  the  ground,  earlier  formation  of  the  planar  blast  wave  leads  to  in¬ 
creased  efficiency.  At  the  limit,  the  distance  between  the  charges  could  be  so  large 
that  DECA  will  not  produce  a  greater  blast  load  than  single  line  charges  of  the 
same  weight  detonated  separately.  That  will  take  place  when  collision  of  the  blast 
waves  from  two  parallel  charges  of  the  array  will  produce  lower  pressure  than  the 
target  threshold.  At  another  extreme,  a  large  amount  of  very  closely  spaced  lines  of 
explosive  will  be  difficult  to  detonate  simultaneously,  and  if  the  line  is  thinner  than 
a  critical  radius  it  will  not  detonate  at  all.  A  feasible  DECA  device  would  have 
parameters  laying  between  these  two  limits,  and  one  of  the  tasks  of  our  proposed 
second-year  development  program  will  be  to  determine  the  range  of  DECA  param¬ 
eters  where  this  concept  is  most  effective.  In  our  current  study,  the  main  task  was 
to  show  the  feasibility  of  DECA,  and  to  correlate  simulations  with  experimental 
tests.  For  this  reason,  we  worked  only  with  two  spacial  distances  between  charges: 
50cm  and  25cm. 

In  order  to  smooth  the  blast  load  distribution  on  the  ground  it  was  first  decided 
to  use  different  explosives.  Iremite  is  a  very  ‘‘slow”  explosive,  with  a  density  of  only 
about  0.9g/cm3.  Its  detonation  velocity  is  3500m/sec.  It  is  also  produced  in  strands 
and  it  could  be  used  is  the  same  way  as  Primacord.  It  was  also  decided  to  use  10 
lines  of  explosive  instead  of  9,  so  the  Z-Y  symmetry  plane  would  lie  between  the 
lines  of  explosive.  That  will  extend  the  size  of  the  array  to  4.5m  x  4m.  In  Figure 
1.6  the  results  of  this  simulation  are  shown  in  the  form  of  pressure  and  impulse 
distribution  on  the  ground.  From  the  graphs  in  Fig.  1.6  it  could  be  concluded  that 


substituting  explosive  from  Primacord  to  Iremite  would  not  improve  the  blast  load 
distribution. 

In  order  to  achieve  a  more  even  blast  load  on  the  ground,  the  distance  between 
was  reduced  from  50cm  to  25cm.  For  the  array  of  4m  x  4m,  25cm  spacing  between 
the  charges  will  lead  to  17  lines  of  explosive.  The  surface  density  of  explosive  in  this 
simulation  was  same  as  in  previouse  simulation  with  Iremite.  In  Figures  l.T(a.b)  and 
1.8(a,b)  contour  plots  of  pressure,  for  X-Z  and  X-Y  cross  sections  corespondingly, 
are  shown  for  t  =  0.15ms  and  t  =  0.35ms.  It  is  easy  to  see  in  these  figures  that  for 
that  case,  the  planar  front  forms  almost  imediatly  after  detonation  and  few  contour 
lines  close  to  the  ground  in  Fig.  1.7a  indicate  that  the  blast  was  originated  from  the 
line  sources.  However,  more  detailed  examination  of  the  pressure  time  history  for 
the  point  under  the  line  explosive,  shown  in  Figure  1.9a,  and  for  the  point  between 
the  lines  of  explosives,  shown  in  Figure  1.9b,  demonstrates  large  variation  of  blast 
load  on  the  ground.  In  Fig.  1.9a  pressure  reached  the  same  maximum  level  as  in 
the  case  of  10  lines  of  Iremite  (see  Figure  1.6),  but  because  the  lines  of  explosive  are 
thinner,  in  the  current  simulation  the  maximum  impulse  is  only  930psi-ms.  Between 
the  lines  the  maximum  pressure  is  4.1kpsi,  twice  as  large  as  in  the  previouse  case, 
and  the  impulse  is  SSOpsi  •  ms. 

Another  method  to  achieve  an  evenly  distributed  blast  load  on  the  ground 
is  to  detonate  DECA  at  some  distance  from  the  ground.  In  that  case,  multiple 
collisions  of  the  blast  waves  from  neiboring  charges  and  formation  of  the  planar 
front  would  begin  prior  to  the  start  of  the  reflection  from  the  ground.  Because 
the  incident  blast  wave  in  this  case  would  be  smoother,  its  reflection  from  the 
ground  will  produce  more  evenly  distributed  load.  In  order  to  prove  this  concept 
we  simulated  detonation  of  the  DECA  aray  of  17  line  charges  of  Iremite  (length 
and  weight  was  the  same  as  in  previous  simulation)  suspended  20cm  above  the 
ground.  In  Figures  1.10(a,b,c,d),  results  for  this  simulation  are  given  in  the  form 
of  pressure,  impulse,  density  and  temperature  time  histories  for  the  points  located 
on  the  ground.  All  these  figures  show  ss  10%  variation  in  pressure  and  less  than 


5%  variation  in  impulse.  That  shows  that  suspention  is  the  most  powerful  practical 
method  to  achieve  even  distribution  of  the  blast  load  on  the  ground.  The  maximum 
pressure  observed  in  simulation  of  the  detonation  of  17  lines  of  Iremite  20cm  above 
the  ground  is  3250psi  and  the  maximum  impulse  was  650psi  •  sec.  This  blast  load  is 
not  suficient  for  the  neutralization  of  an  AT  mine.  In  order  to  increase  the  blast  load 
we  decided  to  use  strands  of  Primacord  explosive,  since  the  Primacord  explosive  is 
about  50 %  more  energetic  than  Iremite. 

In  Figures  1.11  and  1.12  pressure  contours  are  shown  for  the  X-Z  and  X-Y  cross 
sections  of  the  blast  wave  formed  by  detonation  of  an  2.25m  x  2m  DECA  device 
formed  of  10  lines  of  Primacord  spaced  25cm  apart  and  suspended  25cm  above 
the  ground.  Total  weight  of  explosives  in  that  simulation  was  301b.  The  pressure 
contours  in  Fig.  1.11  are  shown  for  t  =  0.06ms,  t  =  0.12ms,  and  t  =  0.27ms. 
Better  grid  resolution  allows  as  to  follow  details  of  the  evolution  of  the  blast  wave 
for  this  DECA  detonation.  At  t  =  0.06ms  the  blast  wave  is  formed  by  cylindrical 
detonations  from  the  parallel  line  charges,  and  colisions  between  the  parallel  waves 
have  only  started.  The  primary  waves  from  the  line  charges  can  be  clearly  located  in 
Fig.  1.11  at  time  O.06ms;  at  that  time  the  front  of  the  blast  wave  has  not  reached  the 
ground.  At  t  =  0.12ms  the  reflection  from  the  ground  begins.  The  reflected  blast 
wave  is  stronger  between  the  charges,  because  of  colision  of  the  blast  waves  from  two 
parallel  sources.  At  t  =  0.27ms  reflection  of  the  blast  wave  is  fully  developed.  As  it 
could  be  seen  in  Fig.  1.12  at  that  time  pressure  is  distributed  evenly  on  the  ground. 
The  levels  of  pressure  and  impulse  could  be  examined  in  more  details  in  Figures 
1.13(a.b.c.d),  where  the  time  histories  of  pressure,  impulse  and  density  are  shown 
for  selected  points  on  the  ground.  The  maximum  pressures  shown  in  Figs. 1.13  are 
in  the  range  of  7kpsi  to  4.5kpsi  and  the  impulse  %  SOOpsi  •  ms.  These  levels  of 
blast  loads  should  be  suficient  to  destroy  an  AT  mine.  Following  this  conclusion, 
in  our  experiments  with  AT  mines  we  used  a  DECA  charge  of  the  same  geometry 
and  surface  density  as  in  the  last  described  simulation,  and  it  was  suspended  25cm 
above  the  ground. 


We  have  studied  the  dependence  of  maximum  pressure  and  impulse,  produced 
by  the  blast  wave  on  the  ground,  on  surface  density  of  Primacord  explosive.  In 
Figure  1.14  the  results  of  this  study  are  shewn  on  the  log(Pressure)  vs.  log( Impulse) 
graph.  DECA  charges  for  this  study  were  about  4m  x  4m  in  size,  suspended  25cm 
above  the  ground.  On  the  same  graph,  vulnerability  limits  are  shown  for  mines  and 
other  military  targets  and  graphs  representing  dependence  of  M-58  line  charge  and 
FAE  on  surface  density  (volume  density  for  FAE)  are  shown.  From  Fig.  1.1 4  it  is 
possible  to  conclude  that  DECA  gives  distinct  advantages  over  M-58  and  FAE.  The 
a  dvantage  of  DECA  over  FAE  is  in  higher  pressure.  According  to  Fig. 1.14,  DECA 
is  more  efficient  than  M-58,  since  it  can  produce  the  same  pressure  on  the  ground 
using  about  fourth  of  M-5S  weight.  Also,  the  impulse  produced  by  DECA  is  larger 
than  M-58’s  impulse. 


2.  EXPERIMENTAL 


All  the  experimental  work  reported  here  was  done  by  SRI  International  of 
Menlo  Park,  California.  In  this  part  of  our  report,  we  will  use  the  data  and  infor¬ 
mation  form  Ref.  3. 

The  main  objectives  of  the  experimental  work  described  in  this  report  were  to 
provide  data  for  correlation  with  the  numerical  simulations,  and  to  perform  tests 
on  inert  mines  that  confirm  the  feasibility  of  this  approach  to  damage  or  destroy 
antitank  mines. 

Approach  and  Test  Procedures 

The  experimental  technique,  to  meet  the  objectives  stated  above,  was  to  deto¬ 
nate  arrays  of  explosive  placed  on  or  above  the  ground  and  to  measure  the  air- 
blast  overpressure  at  several  points  above  and  to  the  side  of  the  array.  Fig¬ 
ure  2.1  shows  a  schematic  layout  of  a  typical  experiment.  Initially,  the  explo¬ 
sive  arrays  consisted  of  ten  4.5m  lines  of  Iremite  spaced  0.5m  apart  and  placed 
on  the  ground  covering  an  area  4.5m  x  4.5m  (Figure  2.2).  Each  line,  was 
made  of  two  or  three  strands  of  1.125  inch  diameter  Iremite.  Each  Iremite 
strand  weighs  about  0.59kg/m  (l.lSlb/ft)  and  the  total  charge  weight  was  from 
52.7kg  (1161b)  to  79.1kg  (173.91b).  The  surface  density  of  the  explosive  arrays  was 
from  2.6kg/m2  (0.531b/ft2)  to  3.9kg/m2  (O.SOlb/ft2).  The  lines  of  Iremite  were  each 
initiated  by  a  single  3.65m  long  (12ft)  strand  of  25gr/ft  Primacord.  These  lead-in 
strands  were  bundled  together  at  a  single  point  and  initiated  by  a  single  9.1m  long 
(30ft)  strand  of  Primacord.  This  initiation  scheme  produces  a  detonation  wave  in 
each  line.  The  detonation  front  from  each  line  sweeps  along  the  array  simultane¬ 
ously.  The  detonation  velocity  of  the  Iremite  is  3500m/s  (11, 500ft /s).  The  entire 
lead-in  and  array  is  detonated  within  3.1ms  of  initiation. 

For  tests  1  through  6,  the  Iremite  was  placed  directly  on  a  prepared  soil  bed. 
The  soil  bed  consisted  of  a  firm  soil  base  that  had  been  leveled.  A  thin  layer  (about 
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5cm)  of  sand  was  placed  over  the  soil  and  packed  in  place  using  a  hand  tamper.  The 
sand  provided  a  smooth  surface  on  which  to  lay  the  Iremite.  In  later  tests  (test  7). 
the  Iremite  array  was  suspended  0.25  m  (10  inches)  above  the  test  bed.  Columns 
of  styrofoam  were  used  at  4  points  along  the  length  of  each  line  to  ensure  a  known 
uniform  distance  from  the  explosive  to  the  ground. 

For  tests  S,  9,  and  10,  the  Iremite  was  replaced  by  bundles  of  400gr/ft  Prima- 
cord.  The  spacing  between  lines  of  explosive  was  reduced  from  0.5m  to  0.25m  and 
the  array  dimensions  were  reduced  to  2.5m  x  2.5m.  However,  the  total  explosive 
surface  density  (2.6kg/m2  )  was  approximately  the  same.  Therefore,  each  line  of 
explosive  consisted  of  7  strands  of  400gr/ft  Primacord.  Each  line  had  a  lineal  den¬ 
sity  of  0.59kg/m  (the  same  lineal  density  of  a  single  strand  of  Iremite).  With  these 
Primacord  lines  on  0.25m  spacings,  the  surface  density  of  the  array  is  2.4kg/m2 
compared  with  2.6kg/m2  for  the  Iremite  array.  •  . 

The  Primacord  arrays  were  suspended  0.25m  above  the  ground  (tests  S  and  9) 
or  0.5m  above  the  ground  (test  10)  in  the  same  manner  as  the  Iremite  array  in  test 
7.  The  Primacord  arrays  were  initiated  in  the  same  way  as  the  Iremite  arrays. 

Instrumentation  for  these  tests  consisted  of  up  to  6  pressure  gauges  mounted 
in  3  steel  poles  (Figure  2.2)  that  were  buried  in  the  ground  in  the  central  area  of 
the  test  bed.  The  poles  were  10.2cm  diameter  (4in)  pipes  approximately  4.Sm  long 
(15.5ft)  long  that  were  buried  1.5m  (5ft)  leaving  3.3m  (10.5ft)  above  the  ground. 
Two  of  the  poles  were  lm  apart  and  third  pole  was  2.25m  to  one  side  (Figure  3). 
Gauge  adaptors  were  inserted  in  each  pole.  Gauges  could  be  placed  1,  2,  or  3m 
above  the  ground  in  each  pole.  For  the  tests  described  here,  only  gauge  positions 
lm  and  2m  above  the  ground  were  used.  For  tests  1,  2.  and  3,  piezoelectric  pressure 
gauges  (PCB  Model  119)  were  used.  Because  of  large  transverse  accelerations  in 
the  steel  poles,  these  gauges  produced  unacceptable  records.  For  tests  4  through 
10.  piezoresistive  pressure  gauges  (Endevco  Model  S530A)  were  used.  These  gauges 
produced  fully  acceptable  records.  The  frequency  response  of  these  gauges  is  about 
90  kHz.  which  means  that  measurements  of  pressure  rise  times  of  about  2 fxsec  to 
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3^5  are  possible.  Figure  2.3  shows  a  layout  of  the  array  for  Tests  4  and  5  showing 
gauge  positions  and  gauge  numbers,  and  the  explosive  positioning. 

Pressure  data  were  recorded  on  a  high  speed  tape  recorder  that  has  a  frequency 
response  of  SO  kHz.  The  analog  tape  was  played  back  through  a  recording  digital 
oscilloscope.  The  digitized  data  were  stored  on  a  floppy  disk  for  later  reduction  by 


a  computer  to  produce  hard  copy  records. 

Inert  antitank  mines  were  exposed  to  the  blast  in  tests  S,  9,  and  10.  These  last 
three  tests  were  designed  to  establish  whether  such  an  explosive  array  could  disable 
an  antitank  mine.  The  inert  mines  consisted  of  a  mine  body  that  was  filled  with 
candle  wax  to  simulate  the  mine  explosive  (Figure  2.4a).  Before  each  test,  an  inert 
booster  and  detonator  were  placed  in  the  pressure  plate  assembly  on  the  top  of  the 
mine  and  the  arming  cap  was  screwed  into  place. 

On  these  tests,  the  explosive  arrays  were  suspended  either  0.25m  (tests  8  and 
9)  or  0.5m  (test  10)  above  the  ground  (Figure  2.4b).  The  antitank  mine  was  placed 
in  the  center  of  the  array.  The  explosion  array  was  positioned  so  that  the  instru¬ 
mentation  pole  containing  gauges  G5  and  G8  was  one  mine  diameter  from  the  edge 
of  the  center  mine  («s  30cm).  The  other  instrumentation  pole  containing  gauges  G4 
and  G7  was  on  the  edge  of  the  small  array.  The  third  pole  was  not  instrumented. 
The  mine  was  placed  on  the  hard  subsoil,  about  5cm  (2in)  below  the  sand  bed 
surface,  and  sand  was  packed  around  the  mine. 

Test  Results 

Table  1  summarizes  the  results  of  the  10  tests  performed  in  this  program. 
As  mentioned  earlier,  tests  1.  2,  and  3  produced  no  usable  data,  although  these 


tests  did  demonstrate  that  the  instrumentation  poles  and  tests  bed  were  able  to 
withstand  repeat  testing.  Tests  4  and  5  demonstrated  the  reproducibility  of  the 
technique.  Figure  2.5  shows  records  from  these  two  tests.  Other  than  differences 
in  the  peak  pressure  of  up  to  35  %,  the  records  are  similar.  The  differences  in  the 
peak  pressure  can  be  attributed  to  the  spacial  gradient  of  the  wave  fronts  as  they 
pass  the  fixed  gauge  position.  Small  changes  in  positioning  of  the  DECA  explosives 
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on  the  ground  could  lead  to  a  change  in  the  effective  gauge  position,  and  that  in 
turn  could  lead  to  a  change  in  pressure  record,  because  of  the  large  spacial  gradient 
of  the  blast  wave  pressure.  Pressure  gauges  G7  and  GS  are  located  symmetrically, 
and  should  produce  similar  pressure  records.  However,  gauges  G7  and  GS  recorded 
substantially  different  levels  of  pressure.  We  will  try  to  explain  this  difference  in 
the  next  part  of  the  report. 

On  test  6,  the  charge  density  was  increased  by  509c  from  2.6kg/m2  to  3.9kg/m2 
The  pressure  measured  at  the  gauge  7  location  lm  above  the  ground  was  3964kPa. 
which  is  an  increase  of  31%  from  the  pressure  measured  in  test  5.  An  increase 
of  50%  in  pressure  should  not  be  expected  in  this  case,  since  pressure  is  linearly 
proportional  to  blast  energy  only  for  a  planar  blast  wave.  In  our  case,  at  the  distance 
of  lm  above  the  ground,  pressure  is  some  complex  function  of  energy  in  between 
the  laws  for  planar  and  cylindrical  blast  waves.  The  increase  in  pressure  is  even 
smaller  for  gauges  G4  and  G5  for  an  unclear  reason.  More  detailed  measurements 
are  needed  in  order  to  study  the  dynamics  of  the  blast  wave  from  DECA.  and  its 
dependence  on  the  energy  and  geometry  of  the.  array. 

In  test  7,  the  nominal  explosive  array  was  suspended  0.25m  above  the 
ground.  The  effect  of  this  configuration  was  to  decrease  peak  pressure  from  about 
3000kPa  (435psi)  to  about  2S60kPa  (415psi)  at  the  lm  point,  and  from  about 
2500kPa  (363psi)  to  about  2240kPa  (325psi)  at  the  2m  point  (Figure  2.6).  The  rise 
time  of  the  pressure  pulse  was  also  increased  when  the  explosive  was  suspended 
above  the  ground. 

For  tests  S.  9,  and  10.  the  Iremite  was  replaced  by  bundels  of  400gr/ft  Prima- 
cord.  The  spacing  between  lines  of  explosive  was  reduced  from  0.5m  to  0.25m  and 
the  array  dimentions  were  reduced  to  2.5m  x  2.25m.  However,  the  average  surface 
density  of  explosive  (2.6kg  /  m2)  was  maintained.  The  Primacord  explosive  (PETX) 
is  about  50%  more  energetic  than  the  Iremite  explosive,  so  for  the  same  surface 
density,  a  larger  pressure  is  expected  from  a  Primacord  array  than  from  an  Iremite 
array.  This  proved  to  be  the  case,  as  evidenced  by  comparing  the  results  of  tests  7 


and  9.  Test  7  used  Iremite  suspended  0.25m  above  the  ground,  whereas  test  9  used 
the  same  areal  density  of  Primacord  suspended  0.25m  above  the  ground.  The  peak 
pressure  at  the  lm  point  was  2S60kPa  (415psi)  for  the  Iremite  array  and  4S60  kPa 
(705  psi)  for  the  Primacord  array,  a  709c  increase  in  pressure.  The  large  increase 
in  pressure  proves  that  for  a  smaller  spacing  of  the  line  charges  the  planar  wave 
will  form  and  that  will  lead  to  higher  overpressures  per  given  weight  of  explosives. 
Also  included  in  tests  9  and  10  were  inert  antitank  mines.  These  mines  included 
a  mine  casing  that  was  filled  with  candle  wax  to  simulate  the  mine  explosive,  a 
pressure  plate  assembly,  and  an  inert  booster  and  detonator.  The  mines  were  about 
30cm  (11. Sin.)  in  diameter  and  about  15cm  (6in.)  high.  They  were  placed  in  the 
center  of  the  Primacord  array,  which  was  25cm  about  the  ground  surface.  The  mines 
were  centered  between  lines  of  primacord  that  were  spaced  25cm  apart.  The  mines 
were  buried  so  that  their  bases  were  on  the  firm  soil  subgrade  about  5cm  (2in.)  be¬ 
low  the  sand  bed  surface.  The  upper  10cm  (4in.)  of  the  mine  was  above  the  ground 
level. 

On  test  9,  the  mine  was  severly  crushed  and  the  pressure  plate,  the  detonator, 
and  the  booster  assemblies  were  completely  blown  off  of  the  mine.  Figure  2.7  shows 
the  damage  to  this  mine.  On  test  10,  the  explosive  array  was  raised  to  50cm  above 
the  ground.  On  this  test,  the  mine  was  crushed  similarly  to  that  in  test  9,  and 
the  upper  portion  of  the  pressure  plate  and  the  detonator  were  blown  away  from 
the  mine.  The  lower  pressure  plate  and  the  rubber  seal  around  the  pressure  plate 
along  with  the  booster  remained  on  the  mine.  This  damage  can  be  seen  in  Figure 
2.Sa.  When  the  mine  was  examined  after  the  test,  the  booster  and  lower  pressure 
plate  assembly  were  easily  removed  from  the  mine.  The  pressure  plate  pieces,  the 
detonator,  and  the  booster  are  shown  in  figure  2.Sb. 
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3.  Comparison  Between  Numerical  Simulations  and  Experiments 

One  of  the  objectives  of  the  current  study  was  to  compare  results  of  numerical 
simulation  of  DECA  to  blast  experiments  with  DECA  charges.  Because  of  bud¬ 
getary  constrains  we  could  compare  only  pressure  history  in  locations  above  the 
DECA  charge,  where  pressure  has  been  measured.  Only  in  6  experiments  of  the  10. 
were  useful  pressure  measurments  made,  and  data  scatter  for  identical  experiments 
is  in  a  30%  range.  In  this  environment,  only  limited  number  of  valid  comparisons 
could  be  made,  and  in  this  section  we  will  discuss  them. 

First,  simulation  corresponding  to  experiments  4  and  5  will  be  compared  with 
experimental  data.  In  Figure  3.1,  measured  in  test  4  and  simulated  pressure  histories 
are  shown  for  the  location  of  gauge  G7.  The  simulated  pressure  closely  repeats  the 
patern  of  the  blast  wave  observed  in  the  experiment.  In  this  case  the  maximum 
pressure  levels  are  predicted  within  10%  of  the  experimental  values.  The  location  of 
first,  second  and  fourth  pressure  peaks  are  predicted  accurately.  However  the  third 
peak  is  missing  in  the  experimental  data.  Study  of  the  pressure  records  at  points 
of  the  computational  domain  adjacent  to  the  point  G7  shows  that  at  a  distance 
of  10cm  from  the  point  G7,  the  third  peak  disapears.  So,  slight  inaccuracies  in 
the  alignment  of  the  explosive  charges  on  the  ground  could  produce  the  observed 
discrepancies  in  the  pressure  time-history  graph. 

Gauges  G7  and  G8  are  symmetrically  located  relative  to  the  vertical  cross  sec¬ 
tion  of  the  DECA  charge,  and  because  we  were  modeling  only  assymetric  part  of  the 
array,  both  of  these  gauges  will  correspond  only  to  one  point  in  the  computational 
domain.  In  Figure  3.2  we  compare  the  pressure  record  of  the  gauge  GS  obtained 
in  test  4  with  numerically  simulated  pressure  for  the  same  point.  The  character  of 
the  graphs  reperesenting  calculation  and  experiment  are  similar;  however,  simula¬ 
tion  predicts  ss  70%  higher  peak  pressure.  The  difference  between  the  simulation 
and  pressure  record  obtained  by  gauge  GS  is  the  same  as  between  gauges  G7  and 
G8,  which  should  show  the  same  pressure,  since  they  located  symmetricalv.  One 
of  the  reasons  for  this  discrepancy  could  be  the  obliqueness  of  the  blast  wave  in 


the  experiment.  The  tests  4  and  5  was  done  with  Iremit  which  has  relativlv  low 
detonation  velocity  of  3500m/sec  At  that  velocity,  the  detonation  wave  will  reach 
the  center  of  the  DECA  array  in  «  0.57ms.  That  is  comparable  to  as  0.42ms  which 
takes  to  the  blast  wave  to  travel  from  the  ground  to  the  location  of  the  gauge  G7  (or 
GS).  That  means  that  at  the  time  the  blast  wave  will  reach  the  pressure  transducer 
it  will  be  oblique  to  horisontal  plane  at  36  deg  angle.  So,  small  deviation  in  the 
orientation  of  gauges  G7  and  G8  can  lead  to  large  differences  in  measured  pressure. 
At  the  same  time  numerical  solution  currently  do  not  simulate  runing  detonation 
wave  (that  was  proposed  in  our  next  years  program)  and  its  assumes  that  array 
detonates  simultaneously.  As  we  will  see  below  in  case  of  2m  x  2m  Primacord  ar¬ 
ray  the  obliqueness  of  the  blast  front  is  smaller  and  the  described  effect  would  be 
negligable.  However,  we  think  that  more  experimental  and  computational  study  is 
needed  to  assess  in  details  formation  of  the  oblique  blast  fronts  and  its  effect  on  the 
ground  blast  load.  That  subject  is  specialy  important  for  design  of  the  initiation 
schemes  for  the  DECA  charges.  Test  5  was  an  exact  repeat  of  test  4  and  results  of 
pressure  measurment  in  this  test  are  very  similar  to  results  of  test  4.  However,  the 
large  difference  in  pressure  measured  in  gauges  G7  and  G8  reported  for  test  4  was 
also  observed  in  test  5.  That  means  that  this  discrepancy  is  systematic. 

In  test  9,  10  lines  of  Primacord  explosive  was  used.  Spacing  between  the  lines 
was  25cm.  In  Figure  3.3  we  compare  the  pressure  history  measured  by  gauge  GS 
in  test  9  and  numerical)'  simulated  results  for  the  same  point.  For  that  point, 
the  results  of  simulation  and  experiment  are  very  close,  although  the  simulation 
predicted  a  stronger  than  measured  rarefraction  wave  behind  the  first  shock  wave. 
In  Figure  3.4,  simulation  and  experiment  are  compared  at  the  point  where  the 
G7  gauge  is  located.  Please  note  that  for  this  experiment,  gauges  GS  and  G7 
are  not  located  in  symmetrical  poisitions.  As  is  shown  in  Fig.  3.4,  gauge  GS  is 
located  75cm  above  the  center  of  the  array  and  gauge  G7  is  located  at  the  same 
height  above  the  array’s  edge.  In  Fig.  3.4,  the  experimentally  measured  maximum 
pressure  is  about  30%  larger  than  simulated.  Also,  simulation  predicts  two  shock 


fronts,  and  the  experiment  measured  only  one  front  for  the  G7  gauge.  The  reasons 
of  these  discrepancies  are  unclear,  and  we  point  again  to  the  inaccuracies  of  laying 
the  Primacord  strands  and  of  placing  the  gauge.  More  experiments  and  calculations 
are  needed  in  order  to  provide  a  better  data  base  for  comparision.  We  believe  that 
the  qualitative  results  of  simulation  for  gauge  G7  are  correct,  since  at  least  two 
shock  fronts  should  be  present  as  a  result  of  colisions  between  the  blast  waves  from 
adjacent  parallel  lines,  equidistant  from  the  point  of  measurment. 

In  Figures  3.5  and  3.6,  results  for  simulation  and  experiment  are  compared  for 
the  points  where  gauges  G4  and  G5  were  located.  Both  gauges  were  located  at  a 
distance  of  1.75m  above  the  ground.  At  that  distance,  as  a  result  of  the  multiple 
colisions  between  the  shock  waves,  the  blast  wave  propagates  with  one  front.  In 
Figures  3.5  and  3.6,  results  of  experimentally  measured  pressure  and  calculation  are 
in  the  range  of  expected  accuracy  for  experiments  and  calculations  (±15%). 


4.  Summary  and  Conclusions 


As  a  result  of  the  one-year,  DARPA  supported  (program  code  5G10/6G10). 
exploratory  effort  to  investigate  the  Distributed  Explosive  Charge  Array  (DECA) 
we  come  to  the  following  main  conclusions: 

a.  Both  experiments  and  computer  simulations  point  to  substantial  advantages  of 
the  DECA  concept  over  known  mine  clearing  methods. 

b.  The  computer  simulations  have  been  found  to  be  in  general  agreement  with  the 
experimental  data,  although  the  accuracy  of  the  experiments  should  be  improved. 
In  particular,  the  ringing  and  heating  of  the  pressure  gauges  should  be  reduced. 
Future  experiments  should  be  better  instrumented,  in  order  to  provide  a  clearer 
picture  of  the  blast  propagation.  Also  a  larger  experimental  data  base  is  needed 
to  clarify  some  inconsistencies  in  pressure  records  (such  as  differences  in  data  from 
gauges  G7  and  GS  in  experiments  4  and  5).  Measurments  on  the  ground  will  add 
very  important  information  about  pecularities  of  the  reflection  of  the  blast  wave 
created  by  DECA. 

c.  Computer  simulation  indicated  that  120  lb  of  DECA  explosive  will  produce 
targeted  impulse  and  overpressure  on  the  ground  (I  %  lpsi  ■  sec,  P  4000psi).  This 
impulse  and  pressure  is  sufficient  for  AT  mine  destruction,  and  will  clear  mines  over 
a  16m2  area.  That  constitutes  an  improvement  by  more  than  a  factor  of  five  over 
the  closest  competetive  device  (M-5S  line  charge). 

d.  The  ability  of  DECA  to  destroy  AT  mines  was  verified  experimentally.  In 
two  experiments  DECA  charges  were  suspended  at  25  cm  and  50  cm  above  the 
ground.  In  both  of  these  experiments  AT  mines  lying  on  the  ground  were  destroyed. 
This  shows  that  relatively  large  variations  in  the  distance  from  the  target  did  not 
affect  substantially  the  ability  of  the  charge  to  destroy  mines.  This  also  indicates 
that  DECA  will  not  need  to  be  deployed  with  high  precision,  which  will  make 
its  production  cost-effective.  However,  the  threshold  at  which  such  mines  can  be 


destoryed  needs  to  be  determined.  The  effects  of  burying  the  mine  below  ground  or 
the  effects  of  reducing  charge  density  (lower  charge  weight,  increased  charge  spacing, 
increased  charge  standoff,  or  less  energetic  explosive)  are  all  parameters  that  need 
to  be  investigated  to  establish  mine  lethality  for  near  surface  explosive  arrays. 
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Figures  and  Table  Captions 


Figure  1.1  Schematic  layout  of  a  typical  simulation. 

Figure  1.2a  Pressure  contours  for  a  9  lines  4m  x  4m  DEC  A  array  test.  X- 
section. 

Figure  1.2b  Pressure  contours  for  a  9  lines  4m  x  4m  DEG  A  array  test.  X- 
section. 

Figure  1.2c  Pressure  contours  for  a  9  lines  4m  x  4m  DECA  array  test.  X- 
section. 

Figure  1.2d  Pressure  contours  for  a  9  lines  4m  x  4m  DECA  array  test.  X- 
section. 

Figure  1.3a  Pressure  contours  for  a  9  lines  4m  x‘4m  DECA  array  test.  X- 
section. 

Figure  1.3b  Pressure  contours  for  a  9  lines  4m  x  4m  DECA  array  test.  X- 
section. 

Figure  1.3c  Pressure  contours  for  a  9  lines  4m  x  4m  DECA  array  test.  X- 
section. 

Figure  1.3d  Pressure  contours  for  a  9  lines  4m  x  4m  DECA  array  test.  X- 
section. 

Figure  1.4a  Impulse  on  the  ground  at  t  =  0.25msec. 

Figure  1.4b  Impulse  on  the  ground  at  t  =  0.50msec. 

Figure  1.4c  Impulse  on  the  ground  at  t  =  1.00msec. 

Figure  1.5a  Time  history  of  pressure,  impulse,  temperature  and  density 
ground  under  the  explosive. 

Figure  1.5b  Time  history  of  pressure,  impulse,  temperature  and  density 
ground  between  lines  of  explosive. 


Z  cross 


Z  cross 


Z  cross 


Z  cross 


Y  cross 


Y  cross. 


Y  cross 


Y  cross 


on  the 


on  the 


Figure  1.6  Maximum  pressure  and  impulse  distribution  on  the  ground  for  10  lines 
(1161b)  of  Iremite. 

Figure  1.7a  Pressure  contours  for  17  lines  4m  x  4m  DECA  array.  X-Z  cross  section. 
Figure  1.7b  Pressure  contours  for  17  lines  4m  x  4m  DECA  array.  X-Z  cross  section. 
Figure  l.Sa  Pressure  contours  for  17  lines  4m  x  4m  DECA  array.  X-Y  cross  section. 
Figure  l.Sb  Pressure  contours  for  17  lines  4m  x  4m  DECA  array.  X-Y  cross  section. 

Figure  1.9a  Time  history  of  pressure,  impulse,  temperature  and  density  on  the 
ground  under  a  line  of  explosive. 

Figure  1.9b  Time  history  of  pressure,  impulse,  temperature  and  density  on  tin 
ground  between  lines  of  explosive. 

Figure  1.10a  Time  history  of  pressure,  impulse,  temperature  and  density  on  the 
ground  for  array  detonation  20cm  above  the  ground.  X=0cm,  Y=0cm,  Z=0cm. 

Figure  1.10b  Time  history  of  pressure,  impulse,  temperature  and  density  on  the 
ground  for  array  detonation  20cm  above  the  ground.  X=10cm,  Y=0cm,  Z=0cm 

Figure  1.10c  Time  history  of  pressure,  impulse,  temperature  and  density  on  the 
ground  for  array  detonation  20cm  above  the  ground.  X=25cm,  Y=0cm,  Z=0cm 

Figure  l.lOd  Time  history  of  pressure,  impulse,  temperature  and  density  on  the 
ground  for  array  detonation  20cm  above  the  ground.  X=40cm.  Y=0cm.  Z=0cm. 

Figure  1.11  Pressure  contours  for  10  lines  (2.25m  x  2.0m)  Primacord  detona¬ 
tion.  X-Z  cross  section. 

Figure  1.12  Pressure  contours  for  10  lines  (2.25m  x  2.0m)  Primacord  detona¬ 
tion.  X-Y  cross  section. 

Figure  1.13a  Time  history  for  pressure,  impulse,  total  pressure  and  density  on  the 
ground.  X=0,  Y=0,  Z=0 

Figure  1.13b  Time  history  for  pressure,  impulse,  total  pressure  and  density  on  the 
ground.  X=5,  Y=0.  Z=0 


Figure  1.13c  Time  history  for  pressure,  impulse,  total  pressure  and  density  on  the 
ground.  X=10,  Y=0,  Z=0 

Figure  1.13d  Time  history  for  pressure,  impulse,  total  pressure  and  density  on  the 
ground.  X='20,  Y=0,  Z=0 

Figure  1.14  Maximum  pressure  and  impulse  as  function  of  surface  (volume  for  FAE ) 
density  of  explosive. 

Figure  2.1  Prespective  of  experimental  configuration. 

Figure  2.2  Explosive  array  and  insrumentation  poles. 

Figure  2.3  Experimental  setup  showing  instrumentation  poles. 

Figure  2.4  Experimental  configuration  for  antitank  mine  tests. 

Figure  2.5  Comparison  of  tests  4  and  5. 

Figure  2.6  Comparison  of  data  to  show  effects  of  suspending  explosive  above  the 
ground. 

Figure  2.7  Damage  to  antitank  mine  in  test  9. 

Figure  2.S  Damage  to  antitank  mine  in  test  10. 

Figure  3.1  Comparison  of  data  from  gauge  G7  (test  4)  and  simulation, 
a  .  Experimental  data 
b  .  Computer  simulation. 

Figure  3.2  Comparison  of  data  from  gauge  GS  (test  4)  and  simulation, 
a  .  Experimental  data 
b  .  Computer  simulation. 

Figure  3.3  Comparison  of  data  from  gauge  GS  (test  9)  and  simulation, 
a  .  Experimental  data 
b  .  Computer  simulation. 

Figure  3.4  Comparison  of  data  from  gauge  G7  (test  9)  and  simulation, 
a  .  Experimental  data 


b  .  Computer  simulation. 

Figure  3.5  Comparison  of  data  from  gauge  G4  (test  0)  and  simulation, 
a  .  Experimental  data 
b  .  Computer  simulation. 

Figure  3.6  Comparison  of  data  from  gauge  Go  (test  9)  and  simulation, 
a  .  Experimental  data 
b  .  Computer  simulation. 
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1.2b  Pressure  contours  for  a  9  lines  4m  x  4m  DEC  A  array  test. 

'C-Z  cross  section.  | 

> 

37 

K 

r$*ki 


% 

m: 


2.5 


-♦  t.»A- 


NRL 
LCP 

DECA  BLAST  ON  THE  GROUND  WEIGHT-] OOLB 
X-Z  PLANE,  Y-0  CM,  TIME  -  0.49 MSEC 


PRESSURE 


0.0  0.5 


1.0  1.5  2.0  2.5  3.0  3.5  4.0 


figS' 

■  f 

Ril 


y.v1 

i  r>. 


X  AXIS 


I*;:*: 1 

§«*< 


I  *1 

pi 

M 


Aril 

M 

r.'l 

W,  if 


P 

[A1 1 

R’tvl 


& 

I 


Figure  1.2c  Pressure  contours  for  a  9  lines  4m  x  4m  DECA  array  test.  X-Z  cross  section. 
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Figure  1.3a  Pressure  contours  for  a  9  lines  4m  x  4m  DEC  A  arrav  test.  X-Y 
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Figure  1.3b  Pressure  contours  for  a  0  lines  4m  x  4m  DECA  array  test.  X-Y  cross  section. 
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Figure  i.3c  Pressure  contours  for  a  9  lines  4m  x  4m  DECA  array  test.  X-Y  cross  section. 
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Figure  1.3d  Pressure  contours  for  a  9  lines  4m  x  4m  DECA  array  test.  X-Y  cross  section. 
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Figure  1.4a  Impulse  on  the  ground  at  t  =  0.25msec. 
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Figure  1.4b  Impulse  on  the  ground  at  t  =  0.50msec. 
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Figure  1.4c  Impulse  on  the  ground  at  t  =  1.00msec. 
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Figure  1.5a  Time  history  of  pressure,  impulse,  temperature  and  density  on  the  ground 
under  the  explosive. 
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Figure  1.5b  Time  history  of  pressure,  impulse,  temperature  and  density  on  the  ground 
between  lines  of  explosive. 
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Figure  1.9a  Time  history  of  pressure,  impulse,  temperature  and  density  on  the  ground 
under  a  line  of  explosive. 
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Figure  1.10a  Time  history  of  pressure,  impulse,  temperature  and  density  on  the  ground 
for  array  detonation  20cm  above  the  ground. 
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Figure  1.10b  Time  history  of  pressure,  impulse,  temperature  and  density  on  the  ground 
for  array  detonation  20cm  above  the  ground. 
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Figure  1.10c  Time  history  of  pressure,  impulse,  temperature  and  density  on  the  ground  tor 
array  detonation  20cm  above  the  ground. 
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Figure  1.11  Pressure  contours  for  10  lines  (2.25m  x  2.0m)  Primncord  detonation. 
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Figure  1.12  Pressure  contours  for  10  lines  (2.25m  x  2.0m)  Primacord  detonation. 
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Figure  1.13b  Time  history  for  pressure,  impulse,  total  pressure  and  densitv  on  the  ground. 
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Figure  2.3  Experimental  setup  showing  instrumentation  poles. 
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Figure  2.4  Experimental  configuration  for  antitank  mine  tests. 
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ABSTRACT 

A  new  three-dimensinal  hydrodynamic  algorithm  has  been  deveoped  utili 
sing  a  virtual  system  for  data  storage.  The  code  called  FAST3D  solves  par 
tial  differencial  equations  aproximated  by  finite  differences.  Data  are 
stored  on  y  -  z  planes,  so  that  only  three  planes  are  in  memory  at  any 
given  time.  The  choice  of  the  algorithm  together  with  the  data  handling 
system  provide  a  highly  accurate  and  efficient  simulation  tool  for  use  on 
vector  computers.  Sample  calculations  are  provided  to  illustrate  the  ver¬ 
satility  of  the  computer  code. 

1.  Introduction 

Three-dimensional  hydrodynamic  algorithms  must  be  accurate  and  highly 
efficient.  Large  amounts  of  memory  coupled  with  large  amounts  of  central 
processer  time  (CPU)  are  required  for  even  moderately  well  resolved  prob¬ 
lems.  FAST2D  addresses  these  points  as  well  as  obtaining  the  required 
storage  by  implementing  a  virtual  system.  Details  of  the  algorithm  and 


4  ► .  at. 


2.  Design  And  Implementation  of  FAST3D 


We  describe  FAST3D,  a  code  which  solves  the  three-dimensional  hydrody¬ 


namic  aquations  for  the  conservation  of  mass,  momentum  and  er.nergy  of  an 


ideal  fluid  using  a  new  leapfrog  FCT  algorithm.  The  equations  advancing 


the  dependent  variables  are 


+  ?-(pv)  -  0; 


3(ov) 


+  7*(pw)  -  -  7(P  -  p  )  -  (p  -  p  )  G(2); 

a  a 


■f  7*(Ev)  -  7 •  ( Pv )  +  S  -  L. 


To  calculate  the  pressure  we  employ  an  equation  of  state  in  the  form 


P  -  P(e,p), 


1  , 

where  e  ■  E  -  -=•  pv^,  In  Eq.  (2),  the  ambient  pressure  P  and  density  p 


are  subtracted  from  corresponding  time-dependent  quantities  in  the  source 


terms  to  avoid  the  "falling  sky"  problem;  G(s),  the  acceleration  due  to 


gravity,  is  generally  taken  to  be  constant  at  its  nominal  sea-level  value. 


S  and  1  denote  energy  sources  and  losses. 


The  solution  is  found  by  numerically  advancing  the  equations  in 


Cartesian  geometry  on  a  rectangular  variably-spaced  mesh.  The  basic 


finite-difference  technique  is  the  fully  three-dimensional  leapfrog  scheme 


with  r.o  time  splitting.  This  scheme  has  very  low  dissipation,  with  a 


linear  amplification  matrix  of  unity  on  a  uniform  grid,  but  is  subject  to  a 


weak  (grid  separation)  instability  and  large  dispersive  errors  at  short 


wavelengths.  These  difficulties  are  overcome  by  using  the  Flux-Corrected- 


Trar.sport  (FCT)  technique,  in  which  the  final  solution  is  made  up  of  the 


weighted  average  of  a  high-order  scheme  (leapfrog)  and  a  low-order  scheme, 


in  this  case  upwind  differencing,  designed  to  maintain  positivity  (nonotc- 


city)  where  it  is  required  physically.  The  FCT  algorithm  is  implemented  in 


the  cede  by  adding  to  the  leapfrog  convective  terms  a  diffusive  flux  pro- 


W5 


ft 

& 

ft 


■ft 

•ft 


BW 


$ 


ional  to  the  absolute  velocity,  plus  a  ssall  velocity-independent  dif- 


ve  term.  If  F  represents  any  of  the  dependent  fluid  variables,  the 
te-difference  scheme  used  to  approximate  the  left-hand  side  of  Eqs. 
is  given  by 

Ilk.  ■  Fuk  *  H  I(Fu>  i-r(Fu)1+1  1  n-.*i  (Fi+r  v 


-  -i_1  (Fi*  )  +  CC?i+1-  2Fi  +  F  +  other  two  directions  (5) 

the  half  of  the  zones  which  are  active  at  any  given  time  step,  and 

Tj.  3  ^  "IF  ui  'u  ^Fi+l“2Fi  +  Fi-l^  +  0Cher  two  directions  (6) 

’  ^  _ 

the  half  which  are  inactive.  Here  C  is  a  constant  used  to  apply  extra 
usica  (beyond  that  needed  to  ensure  positivity)  in  strong-shock  prob- 


Ar.  antidiffusive  step  is  added  to  the  algorithm  in  timesplit  fashion, 
.-i.-.g  most  of  the  added  diffusion.  To  do  this  we  begin  by  using  the 
sported  diffused  momenta  and  density  to  calculate  an  interim  velocity 


ath  direction: 


—  ,  .td,  td 

u1  -  (ou)  /  o  t . 


■  e  calculate  “rav“  antidiffusive  fluxes  for  the  active  zones, 

*  ■  (  —  u  +  r)  -  Ft^) 

i-1/2  1  ix  i-1  }  Ui+1  i  '* 


'1-1/2  ■  !g  "l-l  -  «  <»iU  -  '1U>‘ 


tr.e  inactive  zones, 


1+1/ 2 


<ll“:  *  c>  (F“i  -  F“>' 


4l-l/2  "  “l  +  C)  <?Cl'1  *  Fi-1>-  (,b) 

the  difference  between  neighboring  values  of  the  transported  diffused 


quantities , 


(10) 


A  m  -  F^d 

*1+1/2  i+1  i  ’ 


le  signs  of  the  raw  fluxes, 


Si+l/2 


-  sign($.. 


we  calculate  the  "corrected"  fluxes  according  to  the  Boris-Book  strong 
flux- Uniting  fornula: 


'i+1/ 2 


Si+l/2  nax[0*  ain(  *1+1/2  ,ai+l/2*i+3/2,  ffi+l/2Ai-l/21 .  (12) 


These  fluxes  are  then  used  in  the  antiduf fusion  stage  to  get  the  new  values 
of  the  dependent  variable: 

new  c  r 

5ijk  “  Fijk  +  ^i_i/2  “  ^i+1/2  +  0Cher  CW°  directions*  (13). 

The  flux-correction  part  of  the  solution  leaves  just  enough  of  the  low- 
order  scheae  to  guarantee  acnctonicity  of  the  solutions.  The  driving 
terns  [right-hand  sides  of  Ecs.  (l)-(3)]  can  be  added  in  at  alaost  any 
pcir.t  in  the  calculation.  Currently  this  is  done  following  the  transport 
and  diffusion  in  Eq.  (5). 

The  aain  region  of  the  calculation  eaploys  equally  spaced  zones  to 
nair.tain  high  acuracy.  The  grid,  which  does  not  vary  in  tiae,  may  be 
stretched  ir.  regions  of  little  activity  or  near  the  edges  to  reduce  the 
influence  of  boundary  approxiaatior.s .  Either  reflecting  or  outflow  boun¬ 
daries  nay  be  applied  in  the  transverse  (y  and  z)  directions  and  a  choice 
of  reflecting,  outflow,  or  periodic  boundary  conditions  is  availabe  in  the 


longitudinal  direction.  A  variable  cine  step  is  used  based  on  the  ainirun 
CEL  licit  in  each  direction  over  the  whole  resh,  it  <  ain[dx/(  v  +  c), 
-'•!(  v  -  c),  lz/(  v  +  c)  ] .  The  ncner.tun  equation  contains  a  gravity 


.-gas  equation  of  state  is  also  available  for  strong  shock  cal- 


|WTO^TOKTOtptPnDronorpntfnon^K^^^  — 


•  or  our  applications,  the  shocks  and  other  unsteady  phenomena  being 
simulated  are  initially  localized  near  one  end  of  the  system,  conventionaly 

riken  as  the  origin  of  the  x  axis.  At  early  times  it  is  unnecessary  to  up¬ 
date  variables  far  from  this  region,  so  the  number  of  active  zones  in  this 
d.rection  is  (XX)'  <  (NX) '  .  At  late  times  (NX)'  increases  until  the  en- 


) 


re 


mesh  is  active. 


> 


» 


V 


rnree-dimensional  hydrodynamic  calculations  require  large  amounts  of 
memory,  typically  more  chan  is  physically  available  on  a  given  computer. 
Tr.is  necessitates  the  use  of  auxiliary  storage  such  as  disk.  The  FAST3D 
algorithm  uses  only  the  neighboring  values  of  the  fluid  variables  to  time- 
advance  a  given  value.  The  data  are  stored  on  y-z  planes,  so  that  only 
three  planes  need  be  in  memory  at  any  given  time  during  the  computation 
.?;g.  1).  Furthermore,  by  allowing  I/O  buffers  for  transferring  planes  to 
and  from  the  disk,  computation  can  be  overlapped  with  I/O  to  minimize  over¬ 
all  computer  time.  This  implement at in  of  the  algorithm  requires  two  passes 
: hr curb  the  auxiliary  disk  file  per  time  step:  one  for  the  leapfrog  tic* 
advancement  and  one  for  the  FC7  correction.  The  parallel  architecture  of 
the  Cray  XM?  series  can  be  utilized  here.  Other  passes  through  the  disk 
file  which  do  r.ot  occur  at  every  time  step  may  be  cade  as  needed,  e.g.,  to 
rrccuce  diagnostics  or  restart/dump  files. 

The  leapfrog  FCt  algorithm  involves  enough  artithcetic  to  keep  the 
:rsr.3ccrt  cart  cf  the  FAST2D  code  compute-bound.  However,  due  to  the  ex¬ 
plicit  nature  cf  the  algorithm,  the  computations  on  the  planes  can  be  vrit- 
i=r.  to  :a.<e  advantage  of  vecotr  hardware  that  exists  on  machines  such  as 
i‘5  lray-1.  Vhen  fully  vectorized,  FAST 3D  runs  on  the  Cray-1  at  about 
•  -  -  3  (zone-step).  The  code  has  also  been  run  an  the  Cray-XN.P/12  at  15us  / 
ppr.e  step). 
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3.  Problems  and  Results 

The  code  vas  tested  on  a  variety  of  sinple  idealized  calculations. 

The  FJLenann  (exploding  diaphragm)  problem  vas  used  since  it  embodies  all 
the  gas  dynamic  discontinuities  one  finds  in  shock  physics.  Figure  2  shows 
the  density  as  a  function  of  distance  sons  tine  after  diaphragm  rupture. 
Initially  the  density  and  pressure  increased  across  the  diaphragm  from  left 
to  right  by  factors  of  10  and  10C,  respectively.  An  ideal  gas  equation 
with  a  gamma  of  1.4  vas  used  to  relate  pressure  to  internal  energy.  Exam¬ 
ination  of  Fig.  2  reveals  an  accurate  simulation  of  the  shock  front  moving 
to  the  left,  the  contact  discontinuity  with  only  a  small  amount  of  smear¬ 
ing,  and  the  smoothly  varying  rarefaction  region  propagating  to  the  right. 

A  mere  difficult  problem  is  presented  to  indicate  the  versatility  of 
the  FAST3T  code.  Shock  tubes  with  moderate  shock  strengths  corresponding 
to  over  pressures  of  approximately  three  atmospheres  have  been  used  to  in¬ 
vestigate  the  physics  of  shock  diffraction  in  the  pressence  of  heated  lay¬ 
ers  and  obstacles.  Helium  gas  layers  are  generally  used  to  simulated  the 
higher  temperature  gas  along  the  vail  of  the  tube.  Rigid  obstacles  are 
constructed  and  fixed  downstream  from  the  diaphragm.  When  the  shock  wave 
encounters  the  heated  layer  (helium)  it  begins  to  propagate  faster  and 
creates  an  additional  contact  surface  (discontinuity  in  density  but  not 
pressure).  Furthermore,  the  shock  is  inclined,  moving  upward,  as  one  moves 
away  from  the  shock  cube  vail.  Figure  3  shows  the  density  and  pressure 
contours  in  the  plane  cf  symmetry  (x  -  z  axis).  Note  the  "toeing  out"  of 
the  shock  front  into  the  heated  layer.  When  an  additional  discontinuity 
such  as  ar.  abstacle  is  placed  in  this  flow  field,  a  more  complex  flow  is 


Three-dimensional,  graphics  provide  a  unique  way  to  vie-. 


vr.at  is 


taking  place.  Graphics  programs  nave  been  created  with  the  same 


virtual 


Else,  provide  a  look  loco  the  features  of  interest  in  the  calculation. 
Starting  from  left  to  right,  one  first  sees  a  density  level  of  o  *  5.0  x 
1C~*  which  represents  the  helium  layer.  The  intersection  of  the  three  axis 
is  the  origin  for  the  cartesian  frame  of  reference.  The  shock  wave  is  mov¬ 
ing  from  back  to  the  front.  The  axis  parallel  to  horizontal  is  the  y  axis 


and  the  vertical  axis  is  the  z  axis.  Second  from  left  is  the  density 
level,  0  »  2.0  x  1 0“ 3 ,  which  shows  the  shock  impinging  on  the  rigid  struc¬ 
ture  and  the  diffraction  of  the'  shock  over  the  structure,  respective!}. 


large-scale  scientific  computers  provide  insufficent  central 


aemory  for  even  moderate  simulations.  The  FAST3D  code  combines  a  state-of- 
the-art  numerical  algorithm  with  storage  designed  for  vector  computers. 
Additional  speedups  are  possible  by  taking  a  account  of  the  parallel  fea¬ 
tures  in  the  code.  The  result  is  the  ability  to  model  gas  dynamic  systmes 
which  previously  were  too  complex  and  to  achieve  accuracies  not  otherwise 


.g.  1.  P.epresentation  of  the  data  storage  system  used  in  FAST3D. 

,g.  2.  Density  and  Pressure  vs.  distance  for  the  Rlemann  problem.  The 
trtical  line  in  the  middle  indicates  the  position  of  the  diaphragm. 


Dens i tv  and  Pressure  contours  in 


the  x  -  z  plane  of  symmetry.  The 


shock  is  advancing  from  left  to  right.  Note  the  thermal  layer  along  the 
bottom  in  the  density  figure. 

Fig.  -.  Three-dimensinal  representation  of  the  shock  on  structure  in  a 
heated  layer.  Z  axis  is  vertical  to  page,  y  axis  is  horizontal,  and  x  axis 
moves  out  of  page.  The  shock  is  moving  from  back  to  front.  The  figure  at 


.s  contoured  with  density  -  5  x  104  g/cm^.  The  middi 


e  os  contour- 
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